Switchable resonant coupling of flux qubits 
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We propose a coupling scheme, where two or more flux qubits with different eigenfrequencies 
share Josephson junctions with a coupler loop devoid of its own quantum dynamics. Switchable 
two-qubit coupling is realized by tuning the frequency of the AC magnetic flux through the coupler 
to a combination frequency of two of the qubits. The coupling allows any or all of the qubits to be 
simultaneously at the degeneracy point and can change sign. 



PACS numbers: 74.50.+r, 85.25.Am, 85.25. Cp 

I. INTRODUCTION 

The scheme by Makhlin et alii for coupling supercon- 
ducting qubits was recently followed by several proposals 
for tunable coupling 2 -! 3 -^^^ between superconduct- 
ing qubitsi 1 ^ Broadly speaking, these approaches couple 
qubits through the exchange of virtual excitations in the 
coupler circuit, in which the energy separation Eq be- 
tween the ground and first excited state is much larger 
than the tunnel splitting A in the qubits. The coupling 
strength J is controlled by tuning the energy of the cou- 
pling circuit, via a magnetic field (if the coupler is a 
SQUID), or via a gate voltage (if it is a Cooper pair 
box). 

In principle, the above approaches provide a DC cou- 
pling between qubits, and they allow the realization of en- 
tangling gates only if the difference |A a — Ab\ between the 
tunnel splittings of the corresponding qubits is smaller 
than the coupling energy J. 

In the opposite limit, when |A Q — Af,| ^> \J\, it 
was shownii that the qubit-qubit interaction can be 
controlled by an external variable-frequency magnetic 
field at the combination frequencies, |A a ± Ab\/H. This 
approach^ 1 - (coupling by using a time-dependent mag- 
netic flux, or TDMF) is advantageous due to the reso- 
nant character of the coupling: in experiments it is of- 
ten easier to produce fast and precise frequency shifts of 
the RF control signal, as opposed to changes in the am- 
plitude of the DC signal. The proposal in Ref. [ll] also 
did not require additional, dedicated coupler circuits. Its 
disadvantage was that at least one of the qubits must be 
biased away from the optimal point, which could make 
its operation more difficult and reduce its decoherence 
time. 

Later on, a combination of the TDMF approach with a 
dedicated coupler circuit led to tunable-coupling propos- 
als in Refs. [Tl and [3, where both coupled qubits could 
be simultaneously at their optimal points* 1 ^ 

In this paper we propose an alternative realization of 
the TDMF coupling, which allows to switch the coupling 
on and off, and to change its sign. Our proposal has 



an advantage over both the approaches of Refs. [TlllT^fl3l . 
and the generalization of Ref. [Tj to the Josephson cou- 
pling (described in Section [V] of this paper) , in simulta- 
neously providing: (1) a coupling for arbitrarily biased 
qubits, (2) a higher coupling energy, (3) enhanced pro- 
tection from the flux noise, and (4) the elimination of the 
parasitic first-order DC coupling. 
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FIG. 1: Schematic diagram of two flux qubits and the three- 
junction coupler circuit. 



II. MODEL 

In our proposal, the coupling circuit is a small- 
inductance superconducting loop (L c — > 0) with three 
Josephson junctions (denoted by a, b, c in Fig. [1]). 
The shared junctions a, b ensure a significantly stronger 
qubit-loop co upling than in the case of purely inductive 
(like in Refs. [l2|.13; : ) or galvanic connection^ A control- 
lable DC coupling in a similar device has been recently 
proposed and realized experimentally^ ' 14 : 15 with the cou- 
pling energy Jdc = 1-7 GHz. 

For the sake of simplicity, and without loss of gen- 
erality, we assume that the junctions a and b have the 
same Josephson energy Ej = <&oI c /2tt. The qubit-qubit 
coupling is realized by the (small) junction c, with the 
Josephson energy aEj <C Ej. The coupler circuit (a, b, 
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and c) has a high plasma frequency, u p ~ ^/8EjEc/h, 
so that its energy-level separation is much larger than all 
relevant characteristic energies (J and A a ,b) of the sys- 
tem. The large Josephson energy and large capacitances 
(C » C c « aC) of the coupling junctions ensure that 
Ec/Ej -C 1. This allows us to neglect their degrees of 
freedom and to consider them as passive elements, which 
convert the bias currents I pa and I p b, produced by the 
persistent currents circulating in the attached flux qubits, 
into the phase shift, and therefore the energy shift, of the 
small Josephson junction c. 



III. COUPLING STRENGTH 

Let us now concentrate on the coupler circuit. The 
action of the two flux qubits on it can be represented by 
the bias currents through the junctions a and b. (Due to 
the dominance of the Josephson coupling, we can disre- 
gard the geometric mutual inductances between them.) 
The qubit-qubit interaction energy is obtained by taking 
into account the total potential energy. The latter is the 
free energy of the coupler plus the work performed by 
the qubits a, b on the coupler circuit to keep their persis- 
tent currents constant.— By making use of the quantiza- 
tion condition for the gauge- invariant phase differences, 
(p a — tpt, + <Pc = — 27T$ c /$o- where $ c is the magnetic flux 
through the coupler loop and $o is the magnetic flux 
quantum, the reduced total energy Ut can be written as 

U t = Ut/Ej = - cos ip a - COS Lp b 
-a cos(27r/ c + cp a - (fib) - i pa ip a + i pb <p b , (1) 

where f c = $ c /$o, ipa = I P a/I c , and i pb = I pb /I c - For 
small values of a, i pa , and i po , this potential forms a well 
with a minimum near the point (tp a , <fb) = (0, 0). This is 
the Hamiltonian of a three-junction flux qubrt^£ with bi- 
ased junctions, which can be reduced to the Hamiltonian 
of a perturbed two-dimensional oscillator 
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+ Ejip 2 + + KE,j(ip- -tp*_) 2 
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+ Ej[<p+(ipb — ipa) — <P-(}pa + ipb)]- 



(2) 



Here <p± = ((p a ± <Pb)/2, P± = ihd/dip±, M + = 
2C($ /27t) 2 , M_ = M + (l + 2a), k = 1 + 2acos(27r/ c + 
2tpl), and cp*_ = — asin(27r/ c )/(l + 2acos(2?r/ c )). (see 
Rcf. Il7| for details). The perturbation is Ej[tp+(i p f, — 

'pa 

) - Lp-{i pa + ipb)}- 

From Eq. ©, the normal frequencies of the coupler 
are uj+ = y/2E,j/M+ and u- = ^/2kEjJmZ. Its 
eigenstates, in the lowest order in i pa , i p b, are prod- 
ucts v I'+, m (</J-|-) v E'_ i „(<y9_) of the eigenstates of the normal 
modes. The first-order correction to the ground state en- 
ergy of the coupler is zero, and the coupling energy is 
determined by the second-order correction: 
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Ej(ipa ipb) 2 



(3) 



Thus, it is evident that the coupling is provided by the 
virtual photon exchange between the qubits and the cou- 
pler. 

Separating the term proportional to i pa i P b in the 
second-order correction in Eq. ([3]), we obtain the cou- 
pling energy 



j k> E 



(4) 



Inserting the definition of k into Eq.[4]the coupling energy 
of the three-junction coupler reads 

t _ aEj cos(27r/ c + 2(^1) . . 

1 + 2 a cos(27r/ c + 2<) lpalpb - [ ' 

This expression obviously corresponds to the cr z a z cou- 
pling in the natural basis of qubit states (see, e.g., 
Eqs. (1,3) in Ref . [U) 



H int (t) = J(f c )a z a ^ b . 



(6) 



Obviously, the coupling $Q allows either one or both 
qubits to be in their optimal points. In the experiment 
in Ref. [l5] such interaction was used to realize a tunable 
DC coupling between qubits a and b, by changing the 
coupler bias f c . The strength and sign of the coupling 
depend on the precise value of f c . As mentioned above, 
for time-domain operations it is often easier to manip- 
ulate the frequency of the AC signal / c (t), rather than 
the amplitude of a DC pulse. We will therefore use the 
TDMF approach initially proposed in RefUll. 



IV. EFFECTIVE COUPLING UNDER TDMF 

Let us first consider the effective coupling for an arbi- 
trary J(fc)- Assuming the harmonic flux dependence, 

fc(t) = vq + v x cos fit, 

for the reduced flux in the coupler circuit, and expanding 
J(/ c ) near i>o, we reduce the Hamiltonian of the system 
to 



H(t) = H + H 1 = 

-1^2 A s crJ + (JdcK) + J'M^cos^aX, 



s—a.b 



(7) 



where J'i^o) is the first derivative of the coupling en- 
ergy, taken at vo. In the interaction representation this 
becomes 



H(t) = H + H 1 {t) 



(8) 



with 
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Hi = [JucM + {J'{yo)vi cos fit)] [a z a cos A a t - a v a sin A a t) (a z b cos A b t - a\ sin A b t). 



(9) 



Assuming A a — A5 > and fJ « A a ^ Aj, we see 
that (after averaging over the fast oscillations) only the 
coupling 



survives, where 

Jac = J'{vo)vi- 
The operator in the square brackets, 




(10) 



(11) 



(12) 



is entangling and therefore can be used to construct uni- 
versal quantum ciruits.— 

Our results are somewhat similar to those of Ref. [TH. 
Let us describe the differences. Due to the Joseph- 
son, rather than inductive, coupling, our approach re- 
alizes larger coupling energies, therefore it allows smaller 
values of the parameter a and, correspondingly, is less 
nonlinear. For example, at a = 0.01 and vq = 0.25 
the DC coupling Jdc(0-25) is close to zero and the AC 
coupling (Eq. (fTTjO J'(0.25) = 2TrEj c i a i b is at a max- 
imum (Fig. [5^). Using the experimental value of the 
DC coupling energy for the device shown in Fig. [TJ 
Jo(0) = 1.7 GHz^ we find the AC coupling energy 
J AC = 10- 2 J'(0.25) w 100 MHz (for the reduced mag- 
netic flux amplitude v\ = 10 -2 ). The DC and AC cou- 
plings can be increased by moving to the highly nonlin- 
ear regime with larger a > 10 _1 , but now the points, 
corresponding to "zero" DC coupling and maximal AC 
coupling, do not coincide (Fig. [5b) . 
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FIG. 2: Coupling energy J (solid line) and its first derivative 
J' (dashed line) as a function of the reduced applied magnetic 
flux f c through the coupler loop, for two parameters 01 — 0.01 
(a) and a = 0.1 (b). The J' is scaled by a factor 2ir. When 
the nonlinear response of the coupler increases, the points of 
"zero" DC coupling and maximal AC coupling diverge from 
the common point f c = 0.25. 



V. TIME-DEPENDENT MAGNETIC FLUX 
COUPLING USING AN ITERMEDIARY 
"QUBIT" 

Now let us substitute the smaller junction by three 
Josephson junctions with sizes much smaller than the 
size of the coupling junctions a, b (sec Fig. [3]). This is a 
generalized model of Ref. [l3l . with the inductive coupling 
replaced by the stronger Josephson one. This allows to 
increase A c and decrease the area of the qubits, improv- 
ing the protection of the system against magnetic flux 
noise. Nevertheless we will see that this scheme is at a 
serious disadvantage compared to the coupling of Fig. []] 
because it leads to a strong DC coupling between the 
qubits (i.e., parasitic DC coupling). 



We can now apply the same approach as in Eq. ([3]) . The 
harmonic approximation of ([2]) is now invalid. Instead, 
the coupling energy is determined by the change of the 
ground state energy of the coupling "qubit" c, 



-\y/[2E Jc {<p c 



tt)F 



Here Ej c = $oI pc /2ir and ip c is the phase difference 
across the "qubit" c. Expanding e(27r/ c + 2</>_) at 
27r/ c to second order, we obtain the potential of a two- 
dimensional linear harmonic oscillator with a new value 
of the constant 



k = 1 
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FIG. 3: Schematic diagram of two flux qubits and the five- 
junction ("qubit") coupler circuit. 

where e = e/Ej is the normalized "qubit" energy. Sub- 
stituting the new k in Eq. ((4]) , we arrive at an expression 
for the coupling energy in the simple form 

T _ d 2 s . , 

dipl 

(13) 

where f c = f c — 0.5. The derivative, J' = dJ/df c has a 
maximum at f cm = A c /4/ pc $ : 

JL^JMfcm) i^ 2 )- (14) 

Near this point, the AC coupling energy depends on the 
external magnetic flux only in the second order. This 
formula is equivalent to the expression (25) in Ref. [ID 
provided that we use the standard normalization for 
currents^ i pa = 27rM ac / pa /$ , ipb = 27rM hc / ph /$ and 
neglect the mutual inductance between qubits. It is 
evident from Eq. p4[) that the DC coupling cannot be 
tuned to zero without switching off the AC coupling, 
and it turns out to be much stronger than the lat- 
ter. This is what we refer to as parasitic DC cou- 
pling. Because of large nonlinearity, the AC magnetic 
flux should be much smaller than A c /2/ pc <f>o. Taking, 



e.g., v\ = 10 2 A c /2I pc $ , the AC coupling energy be- 
comes 

JaC = J' ax "I = 10" 2 J(/cm), (15) 

i.e., Jac = 10~ 2 Jdo More importantly, in our proposal 
the DC coupling can be switched off completely, and the 
AC coupling (100 MHz) is five times stronger than in 
Ref. [H. 



VI. CONCLUSIONS 

We propose a feasible switchable coupling between su- 
perconducting flux qubits, controlled by the resonant RF 
signal. Due to the frequency control, it is particularly 
suitable for time-domain operations with flux qubits. 
The coupling energy 100 MHz can be achieved by ap- 
plying a magnetic flux 10 _2 4>o to the coupler with the 
combination frequency 

cj = \A a ±A b \/h. 

The Josephson coupling allows to minimize the area of 
the devices, thus limiting the effects of the flux noise, 
and the coupler thus can act in an almost linear regime, 
which, in particular, suppresses the parasitic DC cou- 
pling. The resulting interaction term also acts as an en- 
tangling gate and enables the realization of a universal 
quantum circuit. 
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